KUNS-1616 



On the Decelerating Shock InstabiUty of Plane-Parallel Slab with 

Finite Thickness 

"Ch : RYOICHI NISHI^ and HIDEYUKI KAMAYA \2 



1 : Department of Physics, Kyoto University, Sakyo-ku, Kyoto 606-8502, JAPAN 
2 : Department of Astronomy, Kyoto University, Sakyo-ku, Kyoto 606-8502, JAPAN 



^ ! email(RN):nishi@tap. scphys.kyoto-u.ac.jp 

^ ! email(HK):kamaya@kusastro. kyoto-u.ac.jp 

T^j- ■ Apr. 16th, 1999 

(N : 

^ : ABSTRACT 

oo ; 
^ . 

■ Dynamical stabihty of the shock compressed layer with finite thickness is 
^ ■ investigated. It is characterized by self-gravity, structure, and shock condition 
0> . at the surfaces of the compressed layer. At one side of the shocked layer, its 

■ surface condition is determined via the ram pressure, while at the other side the 
p thermal pressure supports its structure. When the ram pressure dominates the 
Q ■ thermal pressure, we expect deceleration of the shocked layer. Especially, in 
^ . this paper, we examine how the stratification of the decelerating layer has an 
^ ■ effect on its dynamical stability. Performing the linear perturbation analysis, a 
J> . more general dispersion relation than the previous one obtained by one of the 

^ . authors is derived. It gives us an interesting information about the stability 

. of the decelerating layer; When the deceleration of the layer dominates the 

self-gravity, the so-called DSI (Decelerating Shock Instability) is efficient, while 
the gravitational instability occurs when the self-gravity works better than the 
deceleration. The length-scales of the two instabilities are the order of the 
width of the decelerating layer. Their growth time-scales are the order of the 
free-fall time in the density of the shocked slab. Importantly, they are always 
incompatible. We also consider the evolution effect of the shocked layer. In the 
early stages of its evolution, only DSI occurs. On the contrary, in the late stages, 
it is possible for the shocked layer to be unstable for the DSI (in smaller scale) 
and the gravitational instability (in larger scale), ordinary. The onset-time of the 
gravitational instability is the order of the free-fall time of the external medium. 
After onset, as shown above, growth time-scale becomes the free-fall time of 
the dense slab, which is much shorter than the free-fall time of the external 
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medium. Furthermore, we find there is a stable range of wavenumbers against 
both the DSI and the gravitational instability between respective unstable 
wavenumber ranges. These stable modes suggest the ineffectiveness of DSI for 
the fragmentation of the decelerating slab. Thus, in various cases, the structure 
of the shocked layer should be determined when its stability is discussed. 

Subject headings: hydrodynamics — instabilities — shock waves — 
stars: formation 



1. INTRODUCTION 

The hydrodynamic shock waves are expected to play important roles in many 
astrophysical phenomena, especially in the formation of the structure. The shock waves are 
induced by supernovae (e.g. Tomisaka 1990), expanding Hn regions (e.g. Elmegreen 1989), 
molecular outflows around protostars (e.g. Nakano, Hasegawa, & Norman 1995), stellar 
winds from early type stars (e.g. Stevens, Blondin, & Pollock 1992), and so on. These 
expanding shock waves sweep up the external matter and the swept matter will become the 
materials of the structure at the next stage (e.g. Kimura & Tosa 1988). Previously, it was 
postulated that such mechanism performs in the galaxy formation era (e.g. Ikcuchi 1981). 
In this paper, paying attention to the finite thickness of the shocked layer (cf. Vishniac & 
Ryu 1989), we reexamine its stability by linear perturbation analysis. Especially, we present 
a general dispersion relation including the effect of the self-gravity of a shocked slab with 
finite thickness. 

The stability of decelerating shock waves have been investigated, especially in the 
last two decades. Vishniac (1983) firstly studies the isothermal shock waves expanding 
spherically in a uniform medium with a thin shell picture. Then, he has found that it is 
overstable against linear perturbations and the smaller wavelength perturbation is more 
unstable, although plane shock waves in a steady state were believed to be stable for 
rippling perturbation (Erpenbeck 1962). His instability is attendant on the decelerating 
shock waves (e.g. Nishi 1992). Then, we shall call it DSI (Decelerating Shock Instability). 
Importantly, after his series of papers (Ryu & Vishniac 1987, 1988, 1991 Vishniac & Ryu 
1989) is published, Grun et al. (1991) has reported the existence of instability of shock 
wave in laboratory. After that, Mac Low & Norman (1993) study the nonhnear evolution 
of the blast waves, and confirm Grun et al.. The self-gravity effect on non-linear stage is 
examined numerically by Yoshida & Habe (1992). Dgani, van Buren & Noriega-Crespo 
(1996) also discuss the stability of bow-shocks by means of the numerical simulations, and 
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the transverse acceleration instability semi- analytically. The evolution of the supernova 
remnant is discussed in Chevalier & Blondin (1995). Vishniac (1994) examines the stability 
of a slab bounded on both sides by shocks, which is linearly stable but nonlinearly unstable. 
His analysis is confirmed by numerical simulations (Blondin & Marks 1996). For the 
extended review, see Vishniac (1995). 

There is an important astrophysical phenomenon related to the decelerating shock 
waves. Elmegreen and Lada (1977) propose a scenario on sequential formation of OB 
star subgroups. According to their scenario, an Hn region formed by a group of OB 
stars expands in a molecular cloud and a shock wave sweeps up the gas. Afterwards, the 
matter swept up by the shock becomes gravitationally unstable and fragmentation occurs. 
Elmegreen and Elmegreen (1978) studies the stabihty of the isothermal layer supported by 
the thermal pressure on both sides and showed that the layer is gravitationally unstable 
and the growth rate has the maximum for the perturbation wavelength several times the 
thickness of the layer. The case of the high external pressure to the self-gravity is examined 
by Lubow & Pringle (1993). Welter (1982) calculates the growth rate of the perturbation in 
the layer which is supported by thermal pressure on one side and by balancing ram pressure 
of the shock on the other side. Voit (1988) investigates the incompressible sheet supported 
by unbalanced two thermal pressures and shows the coupling of the gravitational instability 
and the Rayleigh- Taylor instability. The effect of the magnetic field to the stability of the 
self-gravitating slab is precisely examined by Nagai, Inutsuka, & Miyama (1998). 

However, importantly, the shock-compressed layers around Hn regions are supported 
by the thermal pressure on one side and the ram pressure on the other side, and the two 
pressures are not always balanced. In this work, we expect that the instability related to a 
decelerating shock wave plays an important role. Moreover, we must study the instability 
of the shocked self-gravitating layers, since in the scenario of sequential star formation the 
next generation of OB stars is finally formed by self-gravity. Elmegreen (1989) calculates 
the evolution of an isothermal shocked layer and the growth of the perturbations in the 
evolving dense layer with self-gravity. The estimated growth time-scale is ~ 0.25{Gpe)^^^^, 
and is insensitive to the Mach number. Here is the density of the external (preshock) 
medium. On the contrary, Vishniac (1983) estimates a relatively short scale and small mass 
scale for gravitational collapse, which are essentially determined by postshock density. Nishi 
(1992) 's estimate is resemble to Vishniac (1983). The discrepancy originates from the fact 
that they adopt different approximation in their analyses each other. As pointed out by Mac 
Low & Norman (1993), one of the differences is that Vishniac (1983) assumes trans-sonic 
turbulence in the shell, while Elmegreen (1989) assumes supersonic turbulence with Mach 
number close to the shock velocity, and Mac Low & Norman's numerical work, in which 
gravity is neglected, may support Vishniac's assumption. Moreover, very importantly, the 
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structure of the shocked slab are not incorporated in these analysis including gravity. Then, 
as presented just below, we adopt another assumption of the shocked slab to clarify its 
dynamical stability. 

In the current paper, we investigate the instability of a plane-parallel slab with finite 
thickness, which slab is bound by thermal pressure on one side and by ram pressure of 
the shock on the other side. The pressures on the two sides are dissimilar. We assume 
incompressibility for the equation of state of the post shock layer (i.e. in the slab), though 
wc assume compressible fluid at the shock surface. This assumption of incompressibility can 
be suitable, since highly compressed isothermal fluid behaves just like an incompressible 
fluid as shown by Elmegreen and Elmegreen (1978). Moreover, as demonstrated by Mac 
Low and Norman (1993), this approximation may be indicated by the fact that the DSI 
saturates with transonic transverse flows in the shocked slab. By this assumption we can 
obtain, in addition, an analytical solution to the perturbation equations, with which we can 
reveal the precise physics of the instabilities. 

In the following sections, we perform the linear perturbation analysis, to study the 
stability of incompressible slab without self-gravity of the layer in § 2 and with self-gravity 
in § 3. We discuss the DSI and gravitational instabihty of incompressible slab in § 4. We 
also discuss the fragmentation process of the slab in § 5 and summarize our results in the 
final section. 



2. The DSI of incompressible slab 

In order to investigate the stability of the slab of gas swept up by a shock front, we 
analyze an infinite, plane-parallel shock wave advances to the negative direction of z-axis 
with the velocity of in a uniform medium of density of p^. Wc assume incompressibility 
of the gas except at the shock front and the density in the dense slab being p^, and that 
the dense layer is confined from the rear by hot and tenuous gas with a finite pressure Pj. 
As first step, we ignore the self-gravity of the slab to present the dispersion relation of the 
pure DSI (cf. §2 of Nishi 1992). We assume a coordinate that the unperturbed shock front 
is at 2; = 0. At the boundary z = 0, we impose shock boundary conditions: 

[pvi] = 0, (1) 

[P + fml]^0, [pv^v\\]^0 (2) 

where p and P are the density and the pressure, and v± and V\\ are the velocities of the fiuid 
perpendicular and parallel to the shock surface, respectively. The square brackets denote 
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the difference of the enclosed quantity across the shock front. 

The fluid equations for a general first-order perturbation in the slab are 

- V ■ ipsv^) = 0, (4) 

where P is the pressure of the fluid, and the subscripts and 1 refer to unperturbed and 
perturbed quantities, respectively. All perturbed quantities are proportional to e*('^^+'^*), 
then we can rewrite equations (H) and (H) as 

{iuj + {ud^)}vir, - ik (^y^ = 0, (5) 

{luo + {ud,)}v,, - (^j-^ =0, (6) 

ikvix + dzVu = 0, (7) 

where u{= {pe/ Ps)ys) is Vq in the slab and Vi^ and Viz are the x- and z- components of Vi, 
respectively. 

According to Vishniac & Ryu (1989), we get the following three (approximate) solutions 
with two different assumptions. For first set of them, we assume that udz <C iuj- Then, 
we can ignore the terms in parentheses and find a combination of the two independent 
solutions given by 



Pi = P+e^'^ + P-e-'^", 

k Pj. I. k P_ 
Vlx = -e e" 

u Ps UJ Ps 

ik P+ P- 



Vlz = 

U Ps 



■kz 



-kz 



LO Ps 



(9) 
(10) 



To find third solution, we consider the case that the perturbed quantities have large 
gradient and udz terms are important. By inspection the third solution is 



Via 



ps 

Aexp - 



lUO 



u 



(11) 

(12) 
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Auk 



Viz 



UJ 



■ exp 



u 



(13) 



Since we are interested in only the case that perturbation grows, we assume that 
—iuji 3> H/u where if is a thickness of the shocked slab and —iuji means a growth rate of 
the perturbation. 

We shall determine the boundary conditions. Assuming the very high compressibility at 
the shock surface, we can suppose Ps = PeVg. Thus, for the pressure boundary conditions, 
we integrate the equation of hydrostatic equilibrium (Goldreich & Lynden-Bell 1965), and 
find them to the first order; 



Pi(0) = (1 - P)^psr, - 2p,VsViM. (14) 

Pi(iJ) = (l-/5)^p,C, (15) 

where j3 (= Pt/Ps) is the pressure ratio at the surface of the slab, Uo is a surface density 
of the slab, rj{= Zgi — 0) and C(= Zt\ — H) are the displacements of the shock and traihng 
surfaces, respectively. With the shock boundary conditions, we set the following boundary 
condition 

viM = -VsvUO)- (16) 

00 

Applying the boundary conditions, we have 

+ p_e-'^H ^ ^1 _ P)j^kH (P+ - P.e-^'^^) , (17) 

P P f P P \ 

^ + — = (1 - f3)-^kH ^^- — -^Au], (18) 

Ps Ps H^UJ^ \ps Ps J 

where = uVs = peVg/ps. 

Then we find the dispersion relation as 

^k^H^ . „rP/i2 ,1/2 



.2 



C 



2 



± kni^— (1-/9) coth[kH]kH + (1 - [3]^^ 



(20) 



The numerical estimate of (^) is presented in figure 1, where the growth rate is normalized 
by c/H, that is, Q = uj/{c/H). The wave-number is expressed as K = kH and j3 = 0.5 is 
adopted in the same figure. As clearly shown in this figure, we can find that there is a set 
of DSI-mode (the other set is damping mode). 
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We evaluate the limiting case of long wave length {kH — > 0). In the case of — /3) ^ 0, 
we have the growth rate as 

r ^2-'/^(3'/\l-(3y/\kHy/\ (21) 



c/H 



This growth rate is equal to the long wave length limit of equation (2.35) of Nishi(1992) 
with the assumption oi H = ao/ Ps- 



3. The DSI of incompressible slab with self-gravity 

With the effect of self-gravity, the equation of fluid motion for a general first-order 
perturbation is 

where ipi is a perturbed component of the gravitational potential. All perturbed quantities 
are also set to be proportional to e**^'^^+'^*) in this section. Then we can rewrite equations 
with ■^i as 

{iuj + {ud,)}vi:, + ikUi-—\ = Q, (23) 



p. 



s 



{iuj + {ud,)}vu + d,^^r-j^=Q, (24) 

ikvix + dzViz = 0. (25) 

As performed in the previous section, we assume that ud^ <^ ioj firstly. Namely, we ignore 
the advection terms then find a combination of the two independent solutions given by 

V^i = + V-e-'=^ (26) 

Pi = P+e^' + P_c-'=^ (27) 

v.. = f ^ - ij\ e^^ - ^ f ^ - A e-'=^ (28) 
vu ^'^(E±_ ^\ _^1(E^_ e-... (29) 

^ \Ps J ^ \Ps J 

Next, since we are concerned for the large gradient of perturbed quantities, we must 
treat the li^^-terms. By inspection, this third solution is found as 
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Via 



Viz 



Ps 

= Aexp 
Auk 



u 



exp 



u 



(30) 
(31) 
(32) 



Since we are interested in only the case that perturbation grows, we assume that 
—iuji ^ H/u again. 



We must determine the boundary conditions. We have the perturbed column densities 
induced by the displacement of the surface as 



a, = -p,7^e-^"e*('^"+'"*). (33) 
at = p,Ce"^'e*^^"+"*\ (34) 

Since the fluid is incompressible, we can determine potential perturbation, t/^i, from the 
surface density distributions and at alone. With the constraints rj, ( <^2iT/k, we have 
the potentials due to and at as 

= -^!^^e-^^e^(^^+-*), (35) 

rC 

k 

From the equations (PBD, (pSj), and (^) we find 

■At = ^Ce-'«. (37) 
i-- = (38) 



For the pressure boundary conditions, we find, to first order. 



Pi(0) = -2nGplHri + (1 - P)—psV - 2p,VsViz{0), (39) 
Pi{H) = 27:GplHC + (1 - P)—PsC- (40) 
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Applying the shock boundary conditions gives the next boundary condition 



00 



(41) 



Applying the boundary conditions, we get 

27rGps ifP^ 



Ps 



Ps 



.Ps 
V ps 



-2kH 



iAu 



Ps Ps 



p 






(« - 1 


ps 














j 



^ y\ps 

2nGpskH \ (P+ 



Ps 



Ps 



A 



Wsk^ i (P+ 



.Ps 



— iAu 



^ {\Ps ) \Ps 



(42) 

(43) 
(44) 
(45) 
(46) 



where a(= (1 — [3)Ps/ i^iiGaQ)) is the ratio of the deceleration of the whole slab to the 
acceleration of the slab's self-gravity at the surface. 



From equations (^4 ) to (^51) , we have 



(a - 1 
2nGps 



) q {B+-B_-tAu) 

{B_- B_e~^''^ + iAu), 



(47) 



B, + S_e" 



-2kH 



^^(i?+-i?+e-2'=^ + ^Au), 



(48) 



where B^ 
(m) and 



P+l Ps~ V'+ cind i?_ = P^l Ps — ip-i respectively. Then we can rewrite equations 
81) as 



2nGps 



[a 



l)kH^B^ 



+ 



r 00' 



I 2ttGp, 



+ (a-l)A;/7+(l-e-2'=^)|5_ 



+ \l + {a- l)kHjiAu = 0, 



(49) 
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2 

+ e-2'=^(-^ + (a + l)fci7l5_+e-2'=^zv4M = 0, (50) 
and equation (^61) as 

c^k'^B+ - c^A;25_ + (cj^ - c^A;=^)iAM = 0. (51) 
Thus, from equations (^91) ~ (^l]), we have dispersion relation 

uj^ + 27rG'p, |2-^^A;2ff2_2coth[A;i7]A;i/|a;^ 



+ {27iGpsf^{a + l)Y^coth[kH]k^H^ 



1-/3 

-2A;iy+ (l-e-2fc^^) Icu^ = 0. (52) 



We neglect the mode of = 0, because of the assumption of —iujj ^ H/u. Finally, 
then, we obtain the dispersion relation : 



2TxGps 1-/3 2 

± [l^-^^ + coth[kH]kH-lj 

-(1 + a)^ — ^ coth [kH] k^H^ 

~ 1 /2 

+ _^A;2i72 _ (1 _ «2)^2^2 ^ 2A;if + (e-^^^^ - l) | . (53) 



4. Discussions 

First, we investigate the limiting cases of the dispersion relation of equation (^31). To 
the limit of the long wavelength, that is kH — > 0, we have 

- ±kH{l - a^f'\ a = ( a. (54) 



From equation (0), we can classify the characteristics of the instability by a. In the case 
of a > 1, the layer is overstable for long wavelength perturbation. When the self-gravity of 
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the slab does not play important roles, that is, in the limits of a and a oo, we have the 



for gravity modes. 

The numerical presentation of (|53D is given in figures 2-4, where the growth rate is 
normalized by (27rGps)°'^, that is, f2 = uj / . The wave-number is expressed as 
K = kH in all the three figures. In figure 2, we display the cases of a = 1.2 and (3 = 0.2, 
0.4, 0.6 and 0.8. Similarly to figure 1, thus, we find the DSI dispersion relations in figure 
2. For the same value of a, (3 is the larger, the unstable region becomes the narrower and 
the growth rate becomes the smaller. This is because large (5 means inefficiency of the 
deceleration at the surface of the slab. In the case of the adopted parameters, there is not 
instability due to the self-gravity. 

In figure 3, the dispersion relations of a = 0.7 and /3=0.2, 0.4, 0.6 and 0.8 are presented. 
Since a is smaller than unity, we expect the gravitational instability in a range of small 
wavenumbers. Indeed, we find the gravitational instability for all /5 in a range of small K. 
As indicated in this figure, the growth rate of the gravitational instability depends almost 
only on a. On the other hand, we also obtain the DSI in larger K except the case of /3 = 0.8. 
We find again that the growth rate of the DSI does not depend only on a but also on (3. By 
the way, we want to comment here that for all of the three cases {(3 = 0.2, 0.4 and 0.6) is 
stable in the middle of K in figure 3, that is, both of the DSI and gravitational instability 
do not occur for a special range of wave-numbers. For the case of /? = 0.8, moreover, the 
DSI mode is always stable because of small a. These points are the result of the analysis 
incorporating the structure of the shocked slab, then, it does not appear in Nishi (1992). 

In figure 4, we present the dispersion relation for a = a = 0. As clearly shown in the 
solid lines, there is only gravitational instability. In other words, we get no DSI mode. 
Thus, for the case of a, a ^ 0, the self-gravity controls the instability mainly. Indeed, if 
a, a ^ 0, we obtain 



This is the same dispersion relation given by Goldreich and Lynden-Bell (1965). Thus, so 
long as a is different from unity somewhat, the order of the growth rate of the gravitational 
instability becomes [Gps]"^^'^ and the scale-length of the most unstable wavelength becomes 
H (see also figure 3), which are almost in agreement with the result of Elmegreen and 
Elmegreen (1978) 's high pressure case. To obtain this mode, it is important to incorporating 
the structure of the slab. 

Next, we display the linear displacements, rj and C, each of which means rj = Zgi — 



dispersion relation of (pOD. On the other hand, in the case of a < 1, the layer is unstable 




(55) 
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and C = ^ti ~ H, respectively. Figure 5 shows evolution of the linear displacements for the 
case of a ^ 1.0 and K = 0.24. The solid line is Zti and the dotted line is Zgi. From the 
upper to lower panels, we display the evolution of the shocked slab at (a)O.O, (b)l.O ■ 6t, 
and (c)2.0 ■ St where 6t is 7r/4|a;(A; = 0.24if)|. The unit of the length is H. We find the 
bending property of the shocked slab for DSI clearly. On the other hands, figure 6 displays 
the evolution of the gravitational instability of the shocked slab [a, a — > ). Each of the 
displayed time is (a)O.O, (b)l.O ■ 6t, and (c)2.0 ■ 6t, respectively. The one set of evolution is 
very different from that of DSI. Someone might think that the shocked slab has evolved 
nonlinearly. Indeed, the final panels of Fig. 5 and Fig. 6 shows difficulty of the description of 
the linear stage. We present both figures only for reader's clear inspection. The essential 
results of the current discussions are derived from figures 1-4. 

Previously, Yoshida & Habe (1992) discussed the nonlinear stage of DSI and 
gravitational instability under the similar condition adopted in this paper. According 
to them, the gravitational instability has longer wave-length than DSI. This property is 
understood through our linear analysis, since the both modes are incompatible in the 
dispersion relation derived in (0). Thus, our result is in accord with the calculation given 
by Yoshida & Habe. Someone might think that DSI does not have bending shocked layer 
since we can not find it clearly in figure 9 of Yoshida Sz Habe (1992). However, although 
the thickness of our slab is nearly constant in our linear analysis, in the numerical results of 
Yoshida & Habe, the shocked slab has a few times larger thickness than the initial value. 
If the bounding effect were to be more tight in Yoshida & Habe (1992), they would obtain 
our bending property of DSI in their simulations (cf. Chevalier & Theys 1975). 

5. Fragmentation of the shocked slab 

In the previous considerations, we assume that a is always constant in each case. But, 
in reality, a changes with the evolution of the shocked slab (e.g., Elmegreen 1989). In this 
section, then, let's discuss the evolution effects (e.g., Nishi 1992). From the definition, we 
find a (X Ps/ctq (X (Vs/ctq)^. Since the shocked layer sweeps up the matter, the column 
density increases and shock velocity decreases with time, hence a decreases. Thus, the 
character of the instability in the slab changes at the time t = t^r when a becomes unity. 
In the early stages of the evolution of the shocked layer, t < t^r, only the DSI can occur in 
that layer. On the contrary, in the late stages, t > tcv, it is possible for the shocked layer to 
be unstable for the DSI (in smaller scale) and the gravitational instability (in larger scale), 
ordinary. Finally, t ^ tcr, « approaches zero, then the layer becomes stable against the 
DSI and unstable for the gravitational instability only, as equation (55). In the following 
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paragraphs, we shall estimate tcr and discuss evolution of the shocked slab. 

If the shock velocity decreases as Vs (x (—1 < p < 0), we can estimate tcr as follow. 
Column density at the time of t is estimated as 

Vst 

(7o = qPeR where R = — ^ — . (56) 

p + 1 

Here, q is the factor depending on the shape of the shock wave, and, of course, g = 1 is the 
case for the plane shock wave. The distance of the shock surface from the origin is denoted 
as R. From the definition of tcr, (^{t = tcr) = 1, and Ps = peVsi '^^ obtain the following 
relation for tcr'-, 

(p + l)/3V^(l-/3)V^ 

?(27rGpe)V2 ■ ^^'^ 

Thus, tcr-i which is the time for the onset of the gravitational instability, does not depend 
on the initial shock velocity. Although /3 is a function of time, the product of (3^^^ and 
(1 — depends only very weakly on [3. Indeed, for the range of /? = 0.1 — 0.9, which 
is the enough range for realistic situation, fi^^^il — PY^^ has only a value in a range from 
0.55 to 0.71. Thus, even for a realistic situation, the estimation of tcr with the assumption 
of P being constant is fairly valid. That is to say, we estimate tcr as being the order of the 
free-fall time of the external matter. 

In our analysis, the effect of shell expansion is not included. For the case of spherical 
expansion, the shell expansion effect may be important. The onset time of the gravitational 
instability is estimated comparing the divergence rate of the expanding shells to the 
gravitational attraction rate for the expanding shells by Ostriker & Cowie (1981) (in the 
cosmological case) and by McCray & Kafatos (1987) (in the supershell case). Since we 
assume large compression ratio at the shock surface, the estimation for the snowplow phase 
(Eq. (22) of McCray & Kafatos) should be compared with ours. The onset time for this 
case is somewhat shorter than the free-fall time of the external medium for the typical case. 
Thus, it may be correct that tcr is somewhat shorter than the free-fall time of the external 
medium even for the expanding shell case, although more detailed investigation is necessary 
to determine definitely. 

Since the DSI saturates with transonic transverse flows in the shocked slab (Mac 
Low & Norman 1993), only with the effect of the DSI the shocked slab hardly fragments. 
Thus, we should consider the coupling of the gravitational instability with the DSI for 
the fragmentation of the evolving shocked layer. First, we consider the most unstable 
perturbation for the DSI in an early stage {t < tcr)- Considering the constant wavelength, 
the perturbation grows by the DSI, but it may be stable in the later stage as far as the 
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perturbation is linear (see figures 2-4). Q Next, we consider the perturbation being unstable 
for the gravitational instability in the later stages. That perturbation can grow by the 
DSI in the early stages. Then, the DSI might be related to fragmentation process via 
the gravity, which occurs in later stages. By the way, importantly, in the intermediate 
stages, the perturbation becomes stable against both instabilities as shown in figure 3. 
Then, we find that it is difficult for the DSI directly to connect with the gravitational 
instability during the evolution of the slab. The calculation by Yoshida & Habe (1992) 
may suggest the stabilization of the instabilities, that is, the perturbation damps when the 
wavelength is in this stable range probably because of the evolution effect (e.g., Vishniac 
1983, Elmegreen 1989). Thus, as far as the stable range is significant, the effect of the DSI 
for the structure formation may not be important. However, Yoshida & Habe do not get 
the growing bending-mode for all parameters, although Mac Low & Norman clearly present 
the evolution of the bending mode. Thus, there is a possibility that numerical resolution of 
Yoshida & Habe is not enough. Further investigation including nonlinear stage is necessary. 

6. Summary 

In this paper, we present the dispersion relation of the shocked slab with self-gravity. 
According to it, when the deceleration dominates the self-gravity, the so-called DSI is 
efficient. Their growth rates are the function of the scale length which is determined from 
the efficiency of the deceleration. Thus, we confirm previous researches qualitatively. The 
slab becomes gravitationally unstable when the gravitational energy is larger than the 
thermal energy of the layer which is the result of the ram pressure at the shocked surface. In 
other words, when the deceleration is inefficient, the self-gravity works well. Interestingly, 
even if the slab is gravitationally unstable in the long wave-length, it is not gravitationally 
unstable in the short wave-length. Instead of it, the DSI occurs. This suggests that the 
smaller scale structure than the Jeans-length may be result of the DSI. The growth rate 
of the gravitational instability depends on almost only a, while that of the DSI does not 
depend only on a but also on (5. The growth time-scale of the both instabilities is the order 
of the free-fall time for the slab density, while the onset-time of the gravitational instability 
is the order of the free-fall time of the external medium. Thus, within our treatment, the 
fragmentation by the gravitational instability becomes possible only after about the free-fall 
time of the external medium. The analytical new point in this paper is that we find stable 
modes in middle wave-number. This is obtained since we are concerned with the structure 
of the shocked slab with finite thickness (cf. Nishi 1992). This stable range may suggest 



Note that a decreases and iJ increases (i.e. K increases for a constant fc) with time evolution. 
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the ineffectiveness of DSI for the fragmentation of the shocked slab via the evolutional 
effects (Vishniac 1983; Elmegreen 1989; Yoshida & Habe 1992). Thus, wc should determine 
the structure of the shocked layer when we must study the properties of the instabilities 
precisely. 
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FIGURE CAPTIONS 



Fig.l — The numerical expression of The growth rate is normahzed by c/H, that 
is, VL = uj/{c/H). In the figure, fi,. is a real part of fi, while fij is an imaginary part of VL. 
The wave-number is expressed as i^' = kH . The positive real parts are shown in (a), while 



the dotted line corresponds to the plus sign of (^OD and the solid line to its minus sign. In 



(b), the two imaginary parts of both the modes coincide. 



Fig. 2 — The numerical expression of (|53|). The growth rate is normalized by (27rGps)°'^. 
We draw fi,. as a real part of fi, while Vti as an imaginary part of Vt. The wave-number 
is expressed as K = kH. We adopt a = 1.2 and (3 = 0.2(solid), 0.4(dash), 0.6(dot) and 
0.8(dot-dash). Only DSI modes are found for the four examples. 

Fig.3 — The same figure of Fig. 2 except the value of a = 0.7. We also adopt (3 = 
0.2(sohd), 0.4(dash), 0.6(dot) and 0.8(dot-dash). In small K we find the gravitational 
instability, while DSI occurs in large K. In the middle of K, we find just the stable modes. 
The case of (3 = 0.8 does not have DSI even in large range of K. Of course, Qr = for the 
gravitationally unstable modes. 

Fig. 4 — The same figure of Fig. 2 except the value of a, a = 0.0. We find the 
dispersion relation given in Goldreich and Lynden-Bell (1965). The solid line shows the 
gravitational instability, while the dotted line is only for the gravitationally modified sound 
mode. 

Fig. 5 — The one set of the evolution of the linear displacements are displayed. Here, 
■)][= Zsi — 0) and ({= Zfi — H), then we show Zgi as dotted lines and zn as solid lines. Each 
of the displayed time is (a)O.O, (b)l.O ■ 6t, and (c)2.0 ■ 6t. The unit of the length is H. The 
bending property of the shocked layer for DSI is found. 

Fig.6 — The same figure of figure 5. The displayed times are (a)O.O, (b)l.O ■ St, and 
(c)2.0 ■ 6t. The gravitational-unstable slab is presented. 
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